Parallel human genome analysis: microarray-based expression monitoring of 1000 genes.

Microarrays containing 1046 human cDNAs of unknown sequence were printed on glass with high-speed robotics. These 1.0-cm2 DNA "chips" were used to quantitatively monitor differential expression of the cognate human genes using a highly sensitive two-color hybridization assay. Array elements that displayed differential expression patterns under given experimental conditions were characterized by sequencing. The identification of known and novel heat shock and phorbol ester-regulated genes in human T cells demonstrates the sensitivity of the assay. Parallel gene analysis with microarrays provides a rapid and efficient method for large-scale human gene discovery.

PCR-based subtractive hybridization and differences in gene content among strains of Helicobacter pylori

Genes that are characteristic of only certain strains of a bacterial species can be of great biologic interest. Here we describe a PCR-based subtractive hybridization method for efficiently detecting such DNAs and apply it to the gastric pathogen Helicobacter pylori. Eighteen DNAs specific to a monkey-colonizing strain (J166) were obtained by subtractive hybridization against an unrelated strain whose genome has been fully sequenced (26695). Seven J166-specific clones had no DNA sequence match to the 26695 genome, and 11 other clones were mixed, with adjacent patches that did and did not match any sequences in 26695. At the protein level, seven clones had homology to putative DNA restriction-modification enzymes, and two had homology to putative metabolic enzymes. Nine others had no database match with proteins of assigned function. PCR tests of 13 unrelated H. pylori strains by using primers specific for 12 subtracted clones and complementary Southern blot hybridizations indicated that these DNAs are highly polymorphic in the H. pylori population, with each strain yielding a different pattern of gene-specific PCR amplification. The search for polymorphic DNAs, as described here, should help identify previously unknown virulence genes in pathogens and provide new insights into microbial genetic diversity and evolution.

A genomic approach to gene fusion technology

Gene expression profiling provides powerful analyses of transcriptional responses to cellular perturbation. In contrast to DNA array-based methods, reporter gene technology has been underused for this application. Here we describe a genomewide, genome-registered collection of Escherichia coli bioluminescent reporter gene fusions. DNA sequences from plasmid-borne, random fusions of E. coli chromosomal DNA to a Photorhabdus luminescens luxCDABE reporter allowed precise mapping of each fusion. The utility of this collection covering about 30% of the transcriptional units was tested by analyzing individual fusions representative of heat shock, SOS, OxyR, SoxRS, and cya/crp stress-responsive regulons. Each fusion strain responded as anticipated to environmental conditions known to activate the corresponding regulatory circuit. Thus, the collection mirrors E. coli's transcriptional wiring diagram. This genomewide collection of gene fusions provides an independent test of results from other gene expression analyses. Accordingly, a DNA microarray-based analysis of mitomycin C-treated E. coli indicated elevated expression of expected and unanticipated genes. Selected luxCDABE fusions corresponding to these up-regulated genes were used to confirm or contradict the DNA microarray results. The power of partnering gene fusion and DNA microarray technology to discover promoters and define operons was demonstrated when data from both suggested that a cluster of 20 genes encoding production of type I extracellular polysaccharide in E. coli form a single operon.

Genetic analysis using genomic representations

Analysis of the genetic changes in human tumors is often problematical because of the presence of normal stroma and the limited availability of pure tumor DNA. However, large amounts of highly reproducible “representations” of tumor and normal genomes can be made by PCR from nanogram amounts of restriction endonuclease cleaved DNA that has been ligated to oligonucleotide adaptors. We show here that representations are useful for many types of genetic analyses, including measuring relative gene copy number, loss of heterozygosity, and comparative genomic hybridization. Representations may be prepared even from sorted nuclei from fixed and archived tumor biopsies.

Analysis of genomic alterations in benign, atypical, and anaplastic meningiomas: Toward a genetic model of meningioma progression

Nineteen benign [World Health Organization (WHO) grade I; MI], 21 atypical (WHO grade II; MII), and 19 anaplastic (WHO grade III; MIII) sporadic meningiomas were screened for chromosomal imbalances by comparative genomic hybridization (CGH). These data were supplemented by molecular genetic analyses of selected chromosomal regions and genes. With increasing malignancy grade, a marked accumulation of genomic aberrations was observed; i.e., the numbers (mean ± SEM) of total alterations detected per tumor were 2.9 ± 0.7 for MI, 9.2 ± 1.2 for MII, and 13.3 ± 1.9 for MIII. The most frequent alteration detected in MI was loss on 22q (58%). In MII, aberrations most commonly identified were losses on 1p (76%), 22q (71%), 14q (43%), 18q (43%), 10 (38%), and 6q (33%), as well as gains on 20q (48%), 12q (43%), 15q (43%), 1q (33%), 9q (33%), and 17q (33%). In MIII, most of these alterations were found at similar frequencies. However, an increase in losses on 6q (53%), 10 (68%), and 14q (63%) was observed. In addition, 32% of MIII demonstrated loss on 9p. Homozygous deletions in the CDKN2A gene at 9p21 were found in 4 of 16 MIII (25%). Highly amplified DNA sequences were mapped to 12q13–q15 by CGH in 1 MII. Southern blot analysis of this tumor revealed amplification of CDK4 and MDM2. By CGH, DNA sequences from 17q were found to be amplified in 1 MII and 8 MIII, involving 17q23 in all cases. Despite the high frequency of chromosomal aberrations in the MII and MIII investigated, none of these tumors showed mutations in exons 5–8 of the TP53 gene. On the basis of the most common aberrations identified in the various malignancy grades, a model for the genomic alterations associated with meningioma progression is proposed.

Quantitative analysis of chromosomal CGH in human breast tumors associates copy number abnormalities with p53 status and patient survival

We present a general method for rigorously identifying correlations between variations in large-scale molecular profiles and outcomes and apply it to chromosomal comparative genomic hybridization data from a set of 52 breast tumors. We identify two loci where copy number abnormalities are correlated with poor survival outcome (gain at 8q24 and loss at 9q13). We also identify a relationship between abnormalities at two loci and the mutational status of p53. Gain at 8q24 and loss at 5q15-5q21 are linked with mutant p53. The 9q and 5q losses suggest the possibility of gene products involved in breast cancer progression. The analytical techniques are general and also are applicable to the analysis of array-based expression data.

A dynamin GTPase mutation causes a rapid and reversible temperature-inducible locomotion defect in C. elegans

Drosophila shibire and its mammalian homologue dynamin regulate an early step in endocytosis. We identified a Caenorhabditis elegans dynamin gene, dyn-1, based upon hybridization to the Drosophila gene. The dyn-1 RNA transcripts are trans-spliced to the spliced leader 1 and undergo alternative splicing to code for either an 830- or 838-amino acid protein. These dyn-1 proteins are highly similar in amino acid sequence, structure, and size to the Drosophila and mammalian dynamins: they contain an N-terminal GTPase, a pleckstrin homology domain, and a C-terminal proline-rich domain. We isolated a recessive temperature-sensitive dyn-1 mutant containing an alteration within the GTPase domain that becomes uncoordinated when shifted to high temperature and that recovers when returned to lower temperatures, similar to D. shibire mutants. When maintained at higher temperatures, dyn-1 mutants become constipated, egg-laying defective, and produce progeny that die during embryogenesis. Using a dyn-1::lacZ gene fusion, a high level of dynamin expression was observed in motor neurons, intestine, and pharyngeal muscle. Our results suggest that dyn-1 function is required during development and for normal locomotion.

Chromosomal aberrations in PARP−/− mice: Genome stabilization in immortalized cells by reintroduction of poly(ADP-ribose) polymerase cDNA

Depletion of poly(ADP-ribose) polymerase (PARP) increases the frequency of recombination, gene amplification, sister chromatid exchanges, and micronuclei formation in cells exposed to genotoxic agents, implicating PARP in the maintenance of genomic stability. Flow cytometric analysis now has revealed an unstable tetraploid population in immortalized fibroblasts derived from PARP−/− mice. Comparative genomic hybridization detected partial chromosomal gains in 4C5-ter, 5F-ter, and 14A1-C1 in PARP−/−mice and immortalized PARP−/−fibroblasts. Neither the chromosomal gains nor the tetraploid population were apparent in PARP−/− cells stably transfected with PARP cDNA [PARP−/−(+PARP)], indicating negative selection of cells with these genetic aberrations after reintroduction of PARP cDNA. Although the tumor suppressor p53 was not detectable in PARP−/− cells, p53 expression was partially restored in PARP−/− (+PARP) cells. Loss of 14D3-ter that encompasses the tumor suppressor gene Rb-1 in PARP−/− mice was associated with a reduction in retinoblastoma(Rb) expression; increased expression of the oncogene Jun was correlated with a gain in 4C5-ter that harbors this oncogene. These results further implicate PARP in the maintenance of genomic stability and suggest that altered expression of p53, Rb, and Jun, as well as undoubtedly many other proteins may be a result of genomic instability associated with PARP deficiency.

Spectral karyotyping suggests additional subsets of colorectal cancers characterized by pattern of chromosome rearrangement

The abundant chromosome abnormalities in most carcinomas are probably a reflection of genomic instability present in the tumor, so the pattern and variability of chromosome abnormalities will reflect the mechanism of instability combined with the effects of selection. Chromosome rearrangement was investigated in 17 colorectal carcinoma-derived cell lines. Comparative genomic hybridization showed that the chromosome changes were representative of those found in primary tumors. Spectral karyotyping (SKY) showed that translocations were very varied and mostly unbalanced, with no translocation occurring in more than three lines. At least three karyotype patterns could be distinguished. Some lines had few chromosome abnormalities: they all showed microsatellite instability, the replication error (RER)+ phenotype. Most lines had many chromosome abnormalities: at least seven showed a surprisingly consistent pattern, characterized by multiple unbalanced translocations and intermetaphase variation, with chromosome numbers around triploid, 6–16 structural aberrations, and similarities in gains and losses. Almost all of these were RER−, but one, LS411, was RER+. The line HCA7 showed a novel pattern, suggesting a third kind of genomic instability: multiple reciprocal translocations, with little numerical change or variability. This line was also RER+. The coexistence in one tumor of two kinds of genomic instability is to be expected if the underlying defects are selected for in tumor evolution.

Exploring the DNA-binding specificities of zinc fingers with DNA microarrays

A key step in the regulation of networks that control gene expression is the sequence-specific binding of transcription factors to their DNA recognition sites. A more complete understanding of these DNA–protein interactions will permit a more comprehensive and quantitative mapping of the regulatory pathways within cells, as well as a deeper understanding of the potential functions of individual genes regulated by newly identified DNA-binding sites. Here we describe a DNA microarray-based method to characterize sequence-specific DNA recognition by zinc-finger proteins. A phage display library, prepared by randomizing critical amino acid residues in the second of three fingers of the mouse Zif268 domain, provided a rich source of zinc-finger proteins with variant DNA-binding specificities. Microarrays containing all possible 3-bp binding sites for the variable zinc fingers permitted the quantitation of the binding site preferences of the entire library, pools of zinc fingers corresponding to different rounds of selection from this library, as well as individual Zif268 variants that were isolated from the library by using specific DNA sequences. The results demonstrate the feasibility of using DNA microarrays for genome-wide identification of putative transcription factor-binding sites.

Gain of chromosome 3q defines the transition from severe dysplasia to invasive carcinoma of the uterine cervix.

We have chosen tumors of the uterine cervix as a model system to identify chromosomal aberrations that occur during carcinogenesis. A phenotype/genotype correlation was established in defined regions of archived, formalin-fixed, and hematoxylin/eosin-stained tissue sections that were dissected from normal cervical epithelium (n = 3), from mild (n = 4), moderate (n = 6), and severe dysplasias/carcinomas in situ (CIS) (n = 13), and from invasive carcinomas (n = 10) and investigated by comparative genomic hybridization. The same tissues were analyzed for DNA ploidy, proliferative activity, and the presence of human papillomavirus (HPV) sequences. The results show that an increase in proliferative activity and tetraploidization had occurred already in mildly dysplastic lesions. No recurrent chromosomal aberrations were observed in DNA extracted from normal epithelium or from mild and moderate dysplasias, indicating that the tetraploidization precedes the loss or gain of specific chromosomes. A gain of chromosome 3q became visible in one of the severe dysplasias/CIS. Notably, chromosome 3q was overrepresented in 90% of the carcinomas and was also found to have undergone a high-level copy-number increase (amplification). We therefore conclude that the gain of chromosome 3q that occurs in HPV16-infected, aneuploid cells represents a pivotal genetic aberration at the transition from severe dysplasia/CIS to invasive cervical carcinoma.

Exploring drug-induced alterations in gene expression in Mycobacterium tuberculosis by microarray hybridization

Tuberculosis is a chronic infectious disease that is transmitted by cough-propelled droplets that carry the etiologic bacterium, Mycobacterium tuberculosis. Although currently available drugs kill most isolates of M. tuberculosis, strains resistant to each of these have emerged, and multiply resistant strains are increasingly widespread. The growing problem of drug resistance combined with a global incidence of seven million new cases per year underscore the urgent need for new antituberculosis therapies. The recent publication of the complete sequence of the M. tuberculosis genome has made possible, for the first time, a comprehensive genomic approach to the biology of this organism and to the drug discovery process. We used a DNA microarray containing 97% of the ORFs predicted from this sequence to monitor changes in M. tuberculosis gene expression in response to the antituberculous drug isoniazid. Here we show that isoniazid induced several genes that encode proteins physiologically relevant to the drug’s mode of action, including an operonic cluster of five genes encoding type II fatty acid synthase enzymes and fbpC, which encodes trehalose dimycolyl transferase. Other genes, not apparently within directly affected biosynthetic pathways, also were induced. These genes, efpA, fadE23, fadE24, and ahpC, likely mediate processes that are linked to the toxic consequences of the drug. Insights gained from this approach may define new drug targets and suggest new methods for identifying compounds that inhibit those targets.

Knowledge-based analysis of microarray gene expression data by using support vector machines

We introduce a method of functionally classifying genes by using gene expression data from DNA microarray hybridization experiments. The method is based on the theory of support vector machines (SVMs). SVMs are considered a supervised computer learning method because they exploit prior knowledge of gene function to identify unknown genes of similar function from expression data. SVMs avoid several problems associated with unsupervised clustering methods, such as hierarchical clustering and self-organizing maps. SVMs have many mathematical features that make them attractive for gene expression analysis, including their flexibility in choosing a similarity function, sparseness of solution when dealing with large data sets, the ability to handle large feature spaces, and the ability to identify outliers. We test several SVMs that use different similarity metrics, as well as some other supervised learning methods, and find that the SVMs best identify sets of genes with a common function using expression data. Finally, we use SVMs to predict functional roles for uncharacterized yeast ORFs based on their expression data.

Importance of replication in microarray gene expression studies: Statistical methods and evidence from repetitive cDNA hybridizations

We present statistical methods for analyzing replicated cDNA microarray expression data and report the results of a controlled experiment. The study was conducted to investigate inherent variability in gene expression data and the extent to which replication in an experiment produces more consistent and reliable findings. We introduce a statistical model to describe the probability that mRNA is contained in the target sample tissue, converted to probe, and ultimately detected on the slide. We also introduce a method to analyze the combined data from all replicates. Of the 288 genes considered in this controlled experiment, 32 would be expected to produce strong hybridization signals because of the known presence of repetitive sequences within them. Results based on individual replicates, however, show that there are 55, 36, and 58 highly expressed genes in replicates 1, 2, and 3, respectively. On the other hand, an analysis by using the combined data from all 3 replicates reveals that only 2 of the 288 genes are incorrectly classified as expressed. Our experiment shows that any single microarray output is subject to substantial variability. By pooling data from replicates, we can provide a more reliable analysis of gene expression data. Therefore, we conclude that designing experiments with replications will greatly reduce misclassification rates. We recommend that at least three replicates be used in designing experiments by using cDNA microarrays, particularly when gene expression data from single specimens are being analyzed.

Complete genomic sequence of Pasteurella multocida,Pm70

We present here the complete genome sequence of a common avian clone of Pasteurella multocida, Pm70. The genome of Pm70 is a single circular chromosome 2,257,487 base pairs in length and contains 2,014 predicted coding regions, 6 ribosomal RNA operons, and 57 tRNAs. Genome-scale evolutionary analyses based on pairwise comparisons of 1,197 orthologous sequences between P. multocida, Haemophilus influenzae, and Escherichia coli suggest that P. multocida and H. influenzae diverged ≈270 million years ago and the γ subdivision of the proteobacteria radiated about 680 million years ago. Two previously undescribed open reading frames, accounting for ≈1% of the genome, encode large proteins with homology to the virulence-associated filamentous hemagglutinin of Bordetella pertussis. Consistent with the critical role of iron in the survival of many microbial pathogens, in silico and whole-genome microarray analyses identified more than 50 Pm70 genes with a potential role in iron acquisition and metabolism. Overall, the complete genomic sequence and preliminary functional analyses provide a foundation for future research into the mechanisms of pathogenesis and host specificity of this important multispecies pathogen.